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Particle production and correlation functions from Au+Au reactions have been mea- 
sured as a function of both beam energy (2-10.7AGeV) and impact parameter. These 
results are used to probe the dynamics of heavy-ion reactions, confront hadronic models 
over a wide range of conditions and to search for the onset of new phenomena. 
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1. INTRODUCTION 

Very dense nuclear matter is formed in Au+Au collisions at beam energies near 10 
AGeV. The physics of such dense matter is intriguing - whether the system behaves as an 
interacting mixture of excited hadronic resonances, or whether the density is sufficiently 
high that some of the matter is converted to deconfined quarks and gluons. 

Experiments E866 and E917 at the AGS aim to characterize as fully as possible the 
emission of particles from Au+Au collisions to probe the properties of such dense matter. 
In Quark Matter '96[l[j the E866 collaboration reported that the proton rapidity distri- 
bution for central collisions is peaked at mid-rapidity ||. This is consistent with forming 
a dense, baryon-rich system. The dynamics of this system were further characterized by 
the spectra and yield of the most abundantly produced particles, pions and kaons. These 
results are summarized here for reference. The measured mid-rapidity K + /7r + ratio in 
central reactions is 0.19±0.01. This continues a steady increase in the K + //T + ratio from 
p+A@ and Si+A[Q reactions. The transverse m t spectra of different particles have dif- 
ferent spectral shapes, with an overall trend that the proton spectra have larger inverse 
slopes than the kaon spectra, which in turn have larger inverse slopes than the pion spec- 
tra. Detailed analysis [0,0 revealed that the ir~ transverse spectra have a low p t rise above 
an exponential that is steeper than the rise in 7r + spectra. Such features originate from 
the interplay between multiple collisions of hadrons, expansion of the system, decay of 
resonances and Coulomb effects. 

Since Quark Matter '96, the E866 and E917 collaborations have continued to work on 
several open questions, two of which will be the main topics of this paper; 1) is there any 
evidence for new physics beyond hadronic multiple collisions, for example the formation 
of a baryon rich quark-gluon plasma? and 2) is the collision zone at thermal equilibrium? 
These questions are linked because a high collision rate will drive the system towards 
equilibrium and also alter the signatures of a possible quark-gluon plasma that might be 
formed early in the reaction. 

2. EXPERIMENT 

The results in this talk come from two separate collaborations, E866 and E917, that 
used much of the same equipment. Experiment E866 measured Au+Au reactions at 2, 
4, and 10.7 A GeV kinetic energy. Experiment E917 also measured Au+Au reactions 
but at 6, 8 and 10.7 A GeV kinetic energy. Experiment E866 consists of two magnetic 
spectrometers that rotate independently. The Forward spectrometer was used only in 
E866 and was optimized for tracking in the high-multiplicity forward region. This 6msr 
spectrometer consists of a sweeping magnet followed by two tracking stations before and 
after a dipole magnet. Each tracking station is a TPC placed between high resolution 
drift chambers. Particle identification is performed with a TOF-wall with 70ps timing 
resolution. The older spectrometer, Henry Higgins, has an acceptance of 25msr and was 
used by both E866 and E917. It consists of drift chambers on either side of a dipole 
magnet. Particle identification is performed with a TOF-wall with a timing resolution 
of 130ps. The Henry Higgins spectrometer has an on-line particle identification trigger. 
Both experiments used the same detectors for global characterization of the events: a 
multiplicity counter surrounding the target, and a calorimeter placed at zero degrees with 
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a 1.5° opening angle. E917 added a beam-vertexing detector to improve the reaction plane 
measurement. The total systematic uncertainty in the normalization of the cross-sections 
is 15% and the systematic uncertainty on the inverse slope parameters is estimated to 
be 5%. For more details on the equipment and data analysis the reader is referred to 
references [|.|J. 

3. EXCITATION FUNCTION 

Measurement of an excitation function can provide more insight into the mechanism 
of particle production than results at a single beam energy. An excitation function can 
also confront hadronic cascade models over a wider range of conditions and thus permit 
a search for the onset of any new phenomena, e.g. the formation of a small volume of a 
baryon-rich QGP. 
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Figure 1. The invariant yield of K + (left panel) and tt + (right panel) per event as a 
function of m t at mid-rapidity for the different beam energies. The data from E866 are 
shown as filled circles and the data from E917 are shown as open circles. The data at 
10.7AGev are shown at the correct scale, the data at each lower energy are divided by 
ten for clarity. The errors are statistical only. 



Kaon production is a sensitive probe of hadronic multiple collisions in heavy-ion re- 
actions and may also provide evidence for the formation of a QGP. The left panel of 
Figure p] shows a plot of the invariant yield of kaons as a function of transverse mass, 
m t = Jpl + mo, produced in Au+Au central collisions at 2, 4, 6, 8, and 10.7 AGeV ki- 
netic beam energy. The centrality selection at each beam energy was the inner 8% of the 
total interaction cross-section (o"i n t=6.8b). Each spectrum covers the mid-rapidity range 
corresponding to that beam energy, | y ~ Vnn | < 0.25. The kaon spectra were fit with an 
exponential function in m t . The fits reproduce the spectra well and provide the inverse 
slope parameter T and the rapidity density, dN/dy, in that rapidity slice. 
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The right panel of Figure [3] shows the invariant spectra for tt + from Au+Au reactions 
at each of the five beam energies. The pion spectra rise at low p t above an exponential. 
These spectra have therefore been fitted with a double exponential, with the dN/dy and 
mean m t as two of the fit parameters. 

The mid-rapidity yields of tt + and K + as a function of the initial available energy y/s 
are shown in the upper panels of Figure 0. Both pion and kaon yields increase steadily 
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Figure 2. The yield of ir + and K + at mid-rapidity (top-panels) for central Au+Au reac- 
tions as a function of the initial available beam energy. The lower panels show the mean 
m t minus the rest mass for 7r + and K + at the same rapidity. The arrow indicates the N-N 
threshold for K + production. The errors include both statistical and a 5% point-to-point 
systematic uncertainty. 



and smoothly with beam energy. There is no indication of any sudden increase in particle 
yield with increasing y/s. In the lower panels of Figure the mean m t of the transverse 
spectra for pions and kaons are plotted versus y/s. Compared to the increase of particle 
production, the mean m t increases slowly with beam energy. The extra available energy 
mainly goes into more particle production rather than increased transverse energy. 

The increase in kaon yield with beam energy is more rapid than the increase of pion 
yield. This is emphasized in Figure [3], where the K + /tt + ratio is plotted versus y/s. This 
ratio increases steadily from near 3% at 2 AGeV to 19% at 10.7 AGeV. 

These extensive measurements over a broad range of conditions provide a stringent 
test of hadronic models. The predictions of the RQMD hadronic cascade || have been 
compared to the measured ratio K + /n + (Figure^). The model reproduces the data well. 

It is however hard to distill any simple physics lesson from a comparison with such 
a complicated model. One way to improve our understanding is to divide the Au+Au 
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Figure 3. The ratio K + /tt + at mid- 
rapidity in central Au+Au reactions as 
a function of the initial available energy. 
The errors are statistical only. 
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Figure 4. The double ratio K + /tt + at 
mid-rapidity from central Au+Au reac- 
tions divided by K + /tt + of total yields 
from p+p reactions as a function of the 
initial available energy. The errors include 
both statistics and a 5% systematic un- 
certainty in the K/ir ratio from p+p reac- 
tions. 



K + /7r + ratio by the K + /tt + ratio from p+p reactions (Figure [|). This double ratio 
is greater than one, demonstrating that K/ir is enhanced in Au+Au reactions relative to 
p+p collisions. This enhancement is largest at the lowest beam energy, possibly because 
secondary collisions increase in relative importance compared to initial collisions as the 



beam energy is reduced. At the kaon threshold this double ratio goes to infinity JTO] . Also 



note that the double ratio at 2GeV is not plotted because of the large uncertainty in the 
p+p kaon yield. 



4. CENTRALITY DEPENDENCE 

At each beam energy the experiment was triggered primarily with spectrometer-based 
triggers. The spectrometer data can then be sorted off-line into several centrality classes. 
At this stage of the analysis, results are only available for the centrality dependence of 
particle production at the full energy, 11.6 A GeV/c, from Experiment E866. 

As the reactions become more central, the number of primary and secondary hadron- 
hadron collisions increase. By measuring particle spectra and yields as a function of 
centrality we probe how these collisions affect particle production. In addition, because 
the volume of dense nuclear matter increases with centrality JTl]], we can search for the 
onset of any new phenomena associated with high densities. 

The measured energy near zero degrees (E Z cal) is dominated by projectile spectator 
nucleons and provides an estimate of the initial overlap of the projectile and target nuclei. 
The number of projectile participants, N pp , is estimated from the measured Ezcal 

N pp = 197 x (1 - E ZCAL /E^ am ) (1) 
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where E££ m =2123 GeV is the kinetic energy of the beam. The event classes are defined 
by cuts in E Z cal- These are referenced by their cross-section: 0-5%, 5-12%, 12-23%, 
23-39% and 39-73% of the total interaction cross-section, a int =6.8b. 
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Figure 5. The dN/dy and inverse slope distributions as a function of rapidity for K + and 
K~ emitted from mid-central Au+Au reactions at 11.6 A GeV/c. The errors are statistical 
only. The hollow circles are the data points reflected about mid-rapidity (y=1.6). 



In the top panels of Figure [|, dN/dy distributions for K + and K~ are shown for the 
5-12% event class of Au+Au reactions at 11.6 AGeV/c. These have been extracted 
from fits to invariant spectra as a function of transverse mass, similar to those shown 
in Figure [I]. Also extracted from each transverse spectrum is the inverse slope, T. The 
rapidity dependence of the inverse slope parameter is shown in the lower panels of Figure [5[ 
Both the dN/dy and the inverse slope distributions peak at mid-rapidity The dN/dy 
distribution for K~ is narrower than for K + , and the K + inverse slopes tend to be slightly 
larger than the inverse slope for K~. 

To extract the total yield of kaons per event, each dN/dy distribution was fitted with 
a Gaussian centered at mid-rapidity (y=1.6) In Figure ||, the K + yield is plotted as a 
function of the number of projectile participants. The yield of K + increases non-linearly 
with N pp , implying that kaon production is more efficient in central collisions. This is 
emphasized in the middle panel of Figure || where the K + yield per participant is plotted. 
The yield of kaons per participant in central reactions is 3.5 times the yield from isospin- 
averaged N-N collisions at the same beam energy||,[|]. This suggests that the majority of 
kaons in central reactions are produced by secondary collisions. 

The hadronic cascade model, RQMDf6|, attempts to simulate how these secondary col- 
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Figure 6. Left panel: the average yield of kaons in an event as a function of the number 
or projectile participants in Au+Au reactions at 11.6 A GeV/c. Middle panel: the yield 
of K + per projectile participant. Right panel: The ratio of the total yields of K~ and K + 
as a function of the number of projectile participants. The errors are statistical only. The 
hatched boxes are the yields from the isospin averaged initial N-N collisions. 



lisions occur in a heavy-ion reaction. The model reproduces both the absolute magnitude 
of kaon production and the dependence of yield on the number of projectile participants 
quite well. 

The yields of 7r + , 7r~, K~, and p have also been measured as a function of centrality. The 
right panel of Figure |^ shows the ratio of total yields K~ to K + . This ratio is strikingly 
constant over the measured range of centrality and below the estimated K^/K + ratio from 
initial N-N collisions|f||9|]. The K~ yield increases non-linearly with N pp in a manner very 
similar to the K + yield. This is surprising since the two particles have different underlying 
production mechanisms, energy thresholds and absorption rates. 

The observation that the K~/K + ratio is independent of centrality may also test models 
that include an in-medium mass for kaons. There have been several predictions that the 
K~ mass decreases with increasing nuclear density [T^-|T^|. Such a conjecture could lead 
to more phase space for K~+K + pair production which would increase the K~/K + ratio in 
central reactions. This is not observed, potentially placing a constraint on the in-medium 
properties of kaons. 

The role of absorption in heavy-ion reaction and how it affects the final yield of particles 
is perhaps best studied with anti-protons [JL5|]. Figure |7| is a plot of the total yield of p 
per participant as a function of the number of projectile participants. The total yield 
is from a Gaussian fit to the p rapidity distribution measured with the Henry Higgins 
spectrometer [|l6j . 

In contrast to kaon production, the p yield per participant is constant and considerably 
below the N-N value [0]. Both observations are consistent with p being strongly absorbed. 
The final yield of p is sensitive not only to absorption but also to the extent that 
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Figure 7. The yield of p per participant 
as a function of the number of projectile 
participants in Au+Au reactions at 11.6 
A GeV/c. The errors are statistical only. 
The box on the left is the estimated yield 
from the initial N-N collisions. 



Figure 8. The invariant spectrum of p 
at mid-rapidity function of trans- 

verse mass in Au+Au reactions at 11.6 
A GeV/c. For comparison the proton 
spectrum at the same rapidity is also 
shown. The errors are statistical only. 



secondary collisions enhance p production. Extra information to help disentangle these 
two effects may come from the transverse spectra of p. The mid-rapidity spectrum of 



p as measured with the Forward-spectrometer fT6|| in central Au+Au reactions is shown 
in Figure || For comparison the proton spectrum is shown on the same plot arbitrarily 
normalized to the p yield at low m t . Within the current statistics these two spectra have 
a similar shape. This is a puzzle since p produced at low m t should be strongly absorbed, 
which would effectively increase the inverse slope parameter. This effect may, however, 
be counterbalanced by the low phase space available for p production which effectively 
reduces the inverse slope parameter. 

5. PION INTERFEROMETRY 

The interference of identical particles (HBT) provides additional information on the 
dynamics of heavy-ion reactions and on the space-time extent of particle emission. Ex- 
periment E866 has two parallel efforts in this field; 1) to map out the systematic evolution 
of the HBT source parameters with changing initial system size, from peripheral Si+A 
reactions to central Au+Au collisions, and 2) to probe the reaction dynamics by extract- 
ing HBT parameters as a function of the momentum of the particles |IT| , |T8|j . Only the first 
aspect of our efforts will be reported here. 

Experiments E859 (Si+A) and E866 (Au+Au) have measured the systematic evolution 
of HBT interferometry for identical pions. This comprehensive data set establishes how 
the extracted HBT parameters vary as a function of the initial geometry of the reaction. 
Any change in this relationship may indicate the onset of new phenomena, e.g. a change 
in expansion driven by forming a small volume of baryon-rich QGPfl9Hi2T . 
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Each pion two-particle correlation function is measured in three dimensions correspond- 



ing to the three components of relative momentum p9| , |20[| : q z is parallel to the beam, q tout 
is parallel to the sum of the transverse momenta of the pair and qtside is in the orthogonal 
transverse direction. The correlation functions were fit with the following Gaussian. 

C(q) = 1 + \ e ~^ Rtoutqtout }' l + ( R tside(ltside)'' l + {Rzqzf) 1/2 ^) 
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Figure 9. The extracted parameter Rtside (left panel) and the quadratic difference of Rtout 
and Rtside (right panel) versus the cube-root of the total pion dN/dy. E859 Si+Al, Si+Au 
results are shown as open triangles and filled squares respectively. E866 Au+Au results 
are shown as open circles and filled inverted triangles. The errors are statistical only. 



The extracted parameter Rtside is primarily sensitive to the transverse geometry and 
dynamics of the source. In the left panel of Figure || Rtside shows a steady, linear increase 
with the cube-root of the pion rapidity density. There is no anomalous increase that would 
indicate a rapid change in the size of the emitting system or a change in its expansion 
dynamics. 

The quadratic difference of Rtout and Rtside (right panel of Figure ||) is sensitive to both 
the duration of emission of the source and its dynamics |19|J20|| . In contrast to Rtside, the 
quadratic difference is independent of the initial system size. This is consistent with the 
duration of emission being similar for all the systems studied. 



6. TESTS OF EQUILIBRIUM 

It is possible that information from an early, small, volume of QGP may be lost in the 
randomizing collisions during the hadronic phase. Inelastic collisions alter the signatures 
of particle abundances and all collisions alter the kinematic signatures. A system with a 
high collision rate will be driven towards thermal equilibrium and potentially lose infor- 
mation about a possible QGP. To assess this possibility, we need to quantify how close 
the emitting system is to thermal equilibrium. 

As observed earlier in this paper, the rapidity distributions for K + and K~ are different. 
This is emphasized in Figure [K] which shows the K~/K + ratio as function of rapidity for 
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different centrality classes in Au+Au reactions at 11.6 A GeV/c. The rapidity distribution 
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Figure 10. The ratio of the the K~ and K + dN/dy distributions for different centrality 
classes in Au+Au reactions at 11.6 A GeV/c. The errors are statistical only. The lower 
panel is for p+p reactions at 12 GeV/c and is adapted from Fesefeldt et al. M. 



for K~ is narrower than the K + distribution for all centralities and also for p+p reactions 
at 12 GeV/c. This is consistent with K~ having less phase space than K + due to a higher 
energy threshold for K~ production. 

Figure [11] shows the difference between the inverse slopes from K + and K~ transverse 
spectra as a function of rapidity for each centrality class in Au+Au reactions at 11.6 
A GeV/c. Systematically the K + inverse slopes are 10 to 20 MeV larger than the K~ 
inverse slopes. The difference in inverse slopes is consistent with less phase space available 
for K~ production. 

The fact that K + and K~ have different rapidity distributions and transverse spectra 
excludes a commonly used class of thermal models^], namely those models that assume 
that thermal parameters are uniform across the emitting system and that all particles are 
in kinetic equilibrium. For a system at kinetic equilibrium, there is sufficient scattering 
so that the spectra of equal mass particles are identical in shape, even if the system is 
expanding. It may however be possible to reconcile these data with a modified thermal 
model, either with a rapidity-dependent chemical potential or temperature, or possibly 
different freezeout conditions for K~ and K + . 

The comparison between the K + and K~ spectra may also constrain and test hadronic 
models that include in-medium properties of kaons. The effects of a reduced kaon mass 
are often modeled by an attractive kaon mean-field[|r^,0]. This would pull K~ particles to 
lower nif and reduce the inverse slope parameter. The repulsive mean-field for K + should 
increase the inverse slope parameter. Because the mean-fields are predicted to increase 
with density, the difference in kaon spectral shapes might be expected to grow for more 
central reactions. This is not observed. 
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Figure 11. The difference between the K + and K~ inverse slope parameters as a function 
of rapidity for different centrality classes in Au+Au reactions at 11.6 A GeV/c. The errors 
are statistical only. 



7. Summary of Results and Conclusions 

By mapping out the evolution of the spectra, yields and correlation functions of pro- 
duced particles, we have probed the effects of multiple hadronic collisions, absorption, 
and energy thresholds in heavy-ion reactions. Across the broad systematic data set there 
is no evidence for any onset of new behavior beyond hadronic scattering as the beam 
energy or centrality is changed. There are however three outstanding puzzles; 1) why is 
the ratio of yields K~/K + independent of centrality, 2) why do the transverse spectra of 
p and protons have similar shapes, and 3) why is the quadratic difference of Rtout and 
Rtside constant? 

One of the largest puzzles in AGS physics is that hadronic mo dels P, |TT| , p2"| predict that 



a small region of very dense nuclear matter (p > 8po) exists during an Au+Au collision, 
but no evidence for a QGP has been observed. There are two possibilities: 1) a baryon- 
rich QGP is not formed in such a region or is very small, 2) the information from such a 
plasma is lost after hadronization by the many hadronic collisions. The second possibility 
makes it difficult to set a quantitative limit on QGP formation. To do so will require 
a model-dependent study of inserting an ad-hoc change into the particle distributions 
during the model evolution of a heavy-ion reaction. How much information is retained 
from this disturbance provides an estimate of the survivability of QGP signatures. This 
clearly depends on how frequent the collisions are in the hadronic phase. The measured 
spectral differences between K + and K~ can be used to test the assumptions of hadronic 
collision rates in such a study. 

Experimentally, the challenge is to extend the results to more sensitive probes, such 
as the yield of A, <f), and also to extend our studies to more exclusive variables, e.g. the 
azimuthal distribution of particle production with respect to the reaction plane. It will 
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also be important to increase the range of observables included in the excitation function 
to proton stopping via the proton rapidity distributions, as well as to the rapidity and 
m.rp dependence of HBT source parameters. 

This work was supported by the U.S. Department of Energy under contracts with 
BNL, Columbia University, LLNL, MIT, UC Riverside, by NASA, under contract with 
the University of California, by Ministry of Education and KOSEF in Korea, and by the 
Ministry of Education, Science, Sports, and Culture of Japan. 
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